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A d-wave holographic superconductor in the presence of a constant magnetic field is studied 
by perturbation method. We obtain both droplet and triangular vortex lattice solutions. The 
results are the same as that of an s-wave holographic superconductor. The non-Abelian holographic 
superconductor with p + ip-wave background in the presence of a magnetic field is also studied. 
Unlike the d-wave and s-wave models, it is found that the non-Abelian model has only droplet 
£N| ■ solution. 
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I. INTRODUCTION 

The correspondence between a d dimensional quantum field theory and a d+1 dimensional 
gravity theory has provided a new method to understand the strong coupled field theory [l-|4[ . The 
application of this duality to condensed matter physics is helpful to understanding the strong coupled 
many-body systems. Now there are many attempts to use the Gauge/Gra vity correspondence to 
study superfluidity/superconductivity |5l-[l3| (see, for example, Refs. J20l - l23f for reviews). The 
holographic superconductors are possible since there are classic gravity theories in AdS space which 
show local U(l) symmetry broken solutions below a critical temperature [2J, [25[ . Therefore, the 
dual field theories break the global U(l) symmetry and they can be used to study superfluidity or 
superconductivity (in the limit that the U(l) symmetry is gauged). 

The initial holographic superconductor is an s-wave one because the order parameter is a scalar. 
After that follow the non-Abelian holographic superconductors with vector parameters which can 
be dual to a p-wave or a p + ip-wave superconductors. The s-wave holographic model couples the 
Abclian Higgs model to gravity with a negative cosmological constant. One can get solutions which 
spontaneously break the Abclian gauge symmetry via a charged complex scalar condensate near the 
horizon of the black hole when the temperature is low enough. The behavior of an s-wave holographic 
superconductor in the presence of a magnetic field has been studied in many papers [26M33I ] . The 
vortex solution for this model has been constructed in Refs. |29J - [33 |. Especially, Maeda, Natsuume 
and Okamura analytically obtain the same Abrikosov lattice solution as that in the Ginzburg-Landau 



theory [32| . These results indicate that this s-wave holographic superconductor is of type II. The 
coherence length £ is studied in Ref. J28(. £ shows the Ginzburg-Landau behavior (1 — T/T c )~ 1//2 
at the phase transition point. Recently, a <i-wave (spin two) holographic superconductor has been 
constructed, in which the complex scalar field in the s-wave model is replaced by a symmetric 
traceless tensor field whose condensate spontaneously breaks th e ga uge symmetry below T c , and 
becomes zero and so that the symmetry is restored above T c [19l |. We found that the critical 
exponents of the correlation function and the penetration length at T c take the mean-field theory 
values [34||. Another holographic superconductor model of d-wave gap was given in Refs. |35l - [37| . 

The action of the non-Abelian holographic superconductor model consists of SU(2) gauge fields 
and the Einstein-Hilbert action. This Einstein- Yang-Mills (EYM) theory with fewer parameters 
whose Lagrangian is determined by symmetry principles is constructed by Gubser [251 ] and is shown 
to have spontaneous symmetry breaking solutions due to a condensate of non-Abelian gauge fields 
in the theory. Gubser and Pufu studied this model with both p-wave backgrounds and (p + ip)- 
wave backgrounds J17j . Roberts and Hartnoll studied the (p + ip)-wa,ve backgrounds and found 
two major nonconventional features for this holographic superconductor which are different from 
their s-wave counterpart. One is the existence of a pseudogap at zero temperature, and the other 
is the spontaneous breaking of the time reversal symmetry [181 ] . The zero temperature limit of 
the model is studied in Ref. J38|, while in Refs. |39|, E3] the model including back-reactions is 
discussed. In our recent paper |41j , we studied the phase transition properties of this model in the 
presence of a constant external magnetic field. We found that the added background magnetic field 
indeed suppresses the superconductivity. Following closely Maeda and Okamura [28j |. we studied 
the superconducting coherence length and magnetic penetration depth of the p-wave holographic 
superconductor by using perturbation theory near the critical temperature in Ref. |42| . The results 
are the same as the case of the s-wave holographic superconductor which has been studied in Ref. 
|28j . The method that we used to study holographic superconductors is called the "bottom- up" 
approach: we put an arbitrary set of fields into the bulk AdS description and then study the bulk 
theory to get the information on the boundary field theory. It is important to embed the holographic 
superconductors into string theory. In Ref. [43J the authors have studied the top-down approach 
considering various D-brane configurations in the AdS black hole background in the string theory 
framework for the p + ip-wave and s-wave model. Scalars and gauge fields are common in all 
string theory realizations of AdS/CFT and it is reasonable that p + ip-wave and s-wave holographic 
superconductors can be consistently embedded in a string theory setting, whereas the embedding of 
the d-wave model with tensor field into string theory is still an open question [371 ] . 



In this paper, following the method used by Maeda, Natsuume and Okamura in Ref. |32j, we 
analytically study the spatially dependent equations of motion for the <i-wave and p + jp-wave holo- 
graphic superconductor when the added magnetic field is slightly below the upper critical magnetic 
field. The following arc our main results. Firstly, the upper critical magnetic field B c2 for both two 
models is calculated and the phase diagrams are obtained. Secondly, we get the same Abrikosov vor- 
tex lattice solutions for the <i-wave model as that of the s-wave model. Thirdly, for the non-Abelian 
superconductor with p + ip wave backgrounds, we obtain the droplet solutions, but the vortex lattice 
solutions appearing in the s-wave and d-wave models are not possible here. The reason is that the 
Maxwell fields in the non-Abelian holographic superconductor are a subgroup of the 517(2) gauge 
group, hence they do not couple with the condensed fields via covariant derivative like the other two 
models. 

The outline of the paper goes as follows. Section II is devoted to the construction of the triangle 
vortex solution of the d-wave model. In section III we discuss the p + ip-w&ve model's droplet 
solution. Finally, the conclusion and some discussion are given in Section IV. 



II. DROPLET SOLUTION AND VORTEX LATTICE SOLUTION FOR THE d-WAVE 

HOLOGRAPHIC SUPERCONDUCTOR 

In this section we first give the spatial dependent equations of motions for the d-wave model in 
the presence of an uniform magnetic field, then we construct the vortex lattice solution. 

The full gravity theory in 3+1 dimensional spacctime which is dual to a 2+1 dimensional d- wave 
superconductor has the following action [13 ] 






{D^B^YD^B^ + m*B ia ,*B>"' + ^F^F" 



(11.1) 
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where B^ v is a symmetric traceless tensor, R is the Ricci scalar, the 6/L 2 term gives a negative 
cosmological constant and L is the AdS radius, k 2 = 8ttGn is the gravitational coupling. D^ is the 
covariant derivative in the black hole background (D^ = <9 M + iA^ in flat space), q and m 2 are the 
charge and mass squared of B ilVl respectively 

Working in the probe limit in which the matter fields do not back react on the metric as in Ref. 
[lg] and taking the planar Schwarzchild-AdS ansatz, the black hole metric reads (we use mostly plus 
signature for the metric) 



dr 2 



r 



2 



ds 2 = -f(r)dt 2 + f— + —(dx 2 + dy 2 ), (11.2) 



where the metric function f(r) is 



/(r)-^(l-^)- (H.3) 



L and ro are the radius of the AdS spacctime and the horizon radius of the black hole, respectively 
They determine the Hawking temperature of the black hole, 

which is also the temperature of the dual gauge theory living on the boundary of the AdS spacctime. 
Now we introduce a new coordinate z = ro/r. The metric (Eq. (|II.2|l ) then becomes 

dg2 = L 2 a 2 (T) 2 + ^ + dy2 + IW 

z A z z h\z) 



in which h(z) = 1 - z 3 and a(T) = r /L 2 = 4ttT/3. 

For the d-wave backgrounds, the spatial dependent ansatz takes the following form |28| 

B^ v = diagonal (0, 0, f(z, x, y), -f(z, x, y)) , A= cf>(z, x, y)dt + A y (z, x, y). (II. 6) 

We assume the vector potential A y is nonvanishing since we need a non- vanishing magnetic field on 
the boundary. A x can be set to zero when we take a suitable gauge. This ansatz for the tensor field 
captures the feature of a d-wave superconductor in which there is a condensate on the x-y plane on 
the boundary with translational invariancc, and the rotational symmetry is broken down to Z(2) 
with the condensate changing its sign under a ir/2 rotation on the x-y plane. 
With this ansatz, we can derive the equations of motion: 



z or or u or or z a z h z l or z z 



a m\A, + aB.hB.A, + l -SiA, + H!ijCM _ IjiJML _ ** 2 -yj f = o. ,11.9) 

In order to solve the above equations, we have to introduce the following boundary conditions on 
the horizon and the boundary: 

(i) On the horizon(z = 1), the scalar potential <j> = since the 4>dt must be well defined. The 
other fields should be regular. 

(ii) On the boundary (z = ), w e are only interested in the L 2 m 2 = —1/4 case. The boundary 
conditions for /, <fi and A y are [281 ] : 

/ = A z A - + A lZ A + + ■ ■ ■ , (11.10) 

4> = n-pz + --- , (n.ii) 

B(x) = x A y - d y A x , (11.12) 



in which A ± = - 1= F^i7+-L 2 » ii and L 2 m 2 > _^ In the case we stu died A_ = -5/2 and A + = 3/2. 
Ao is the source, then A\ is the vacuum expectation value (VEV) of the operator that couples to 
B in the boundary theory. Aq can be set to zero [19( . The order parameter of the boundary theory 
can be read off from the asymptotic behavior of tensor field B^ y , 

(C V ) - ( t \ ) (IL13) 

where (i,j) are the indexes in the boundary coordinates (x,y). [i is the chemical potential and p is 
the charge density of the field theory. -B(x) is the magnetic field of the field theory on the boundary. 
To exactly solve the above nonlinear coupled partial differential equations is a difficult task. But 
we can perturbatively solve these equations when the magnetic field is slightly below the upper 
critical field B C 2- First we define a small parameter e = (B C 2 — B)/B C 2, then we can expand the 
fields as : 

f(x, z) - e 1/2 /i(^, z) + e 3/2 / 2 0z, «) + •■-, (H.14a) 

A y (x,z)=A^(x,z) + eA^(x,z) + --- , (II.14b) 

<j>(x, z) = (o) (cc, z) + e<f) {l \x, z) + ■ ■ ■ (II.14c) 



in which x = (x,y). The zeroth order solution corresponding to the normal state is 



/ = ,<f, = n(l-z), A°=B c2 x. 



(11.15) 



We can see clearly that the magnetic field on the boundary is B c2 . Substituting Eq. (|II.15[> into 
the equations of motion, with the following masatz fi(x, z) = e tpy F(x, z;p)/L (p is a constant), the 
equation of motion for F is 



hd 2 z + {d z h + — )d z + 

1 



2d z h p 2 {\-zf 



o- 



dx 2 



( P - B c2X y 



a 2 h 
F(x,u;p). 



4ft 



L 2 m 2 f 



F(x,u;p) 



(11.16) 



Then we separate the F as F n (x,z;p) = p n (z)~f n (x;p)/L, where A„ is a constant. p n and j n admit 
the following equations: 



ifl 



dx 2 



x*_ 

4 



ln(x;p) 



A, 



ln{x;p), 



2h s 



hdi + {d z h + —)d zPn (z) 



2r2 



m 2 L 



-{l-zf + 



4h 2d z h 



(II.17a) 



(II.17b) 



where X := y / 2B c2 (x — p/B C 2) q '■= n/& are dimensionlcss, L 2 m 2 = —1/4 in our calculation. Eq. 
(II. 17a) determines the distribution of the order parameter on the x — y plane, while Eq. (II. 17b) 
determines when a superconducting phase transition will occur. 

The solution of (|II.17a|) that satisfies the boundary condition and limui^oo |7 n | < oo can be 
expressed in terms of the Hermitc functions H n as follows 



-X 2 /4 



ln{x]p) 

and the corresponding eigenvalue A„ is 

A„ = 2n + 1 



Hn(X), 



(11.18) 



(11.19) 



where n = 0, 1, 2, 3 • • ■ . The n = solution is the droplet solution, and the vortex solution can be 
constructed from the droplet solution: 



7o(cc;p) = e x ' 4 = cxp 



1 



-P r o) 



(11.20) 



where ro := l/y/B C 2- 

Before we construct the vortex lattice solution, let us discuss the phase diagram. From Eq. 
(II. 17b) we can obtain the phase diagram. The upper critical magnetic field given by this equation 
has a non-zero solution satisfying the boundary conditions. This can be done numerically for a given 
q, which corresponds to a fixed temperature. We can find a critical value B c2 / [M 2 above which the 
equation have only vanishing solution. The maximum upper critical magnetic field is given for n = 
when (A ra ) take the minimum value. In Fig. 1 we give the phase diagram, from which we can find 
that B C 2 ex (1 — B / B c2 ) around T c . This is the same as the BCS theory. 

Since A„ is independent of p, a linear superposition of the solutions e zpv po(u)^/o(x; p) with different 
p is also a solution of the equation of motion for /i : 



/i(*,u) = ^5>e*»»<ybfo Pi)- 

i 



(11.21) 



B c2 /p 2 




FIG. 1. The phase diagram of a a!- wave holographic superconductor in the presence of a magnetic field. 



Here we get the most important result in this section. When we choose a suitable configuration of q 
and pi, we can construct triangular lattice solutions. It is very interesting that the result Eq. (11.21) 
is very similar to the expression of the order parameter of G-L theory for the type II superconductor 
in the presence of a magnetic field when B = B C 2 , which is 



01 



J2cie mv cxp[ 



Xl. 



2e 



(11.22) 



where £ is the superconducting coherence length, 
Eq. (11.21) with Eq. (11.22), we get 



Xi 



A* 



-^ , and $o is the flux quantum. Comparing 



B c2 oc 



r 



(11.23) 



which is also similar to the result of the GL theory. According to the behavior that B C 2 ex (1 — T/T c ) 
near T c , we have £ oc (1 — T/T c ) _1 / 2 . This result is also the same as that of the GL theory. We have 
also obtained this result by another way in Ref. [34| . 

Thus, the construction of triangular lattice from droplet solutions is similar to what Abrikosov 
did in his initial paper. This procedure has been made for the s-wave model in Ref. [32|]. In the 
d-wave model, the construction process is the same. We briefly review the result below, considering 
the following form of pi and c; : 






L 



l — ~oo 



Cl 



exp 



K< 2 ) 



pi 



2irl 
air ' 



(II.24a) 
(II.24b) 



for arbitrary parameters a\. The solution in Eq. (|II.24[) represents a lattice. o(x) := |7i(a;)| 2 in 
which the fundamental region Vq is spanned by two vectors bi — air d y and fo 2 = 27rro/ai9 x + 
airo/2d y , and the area is given by 27T7-Q. Then the magnetic flux penetrating the unit cell is given 
by B C 2 x (Area) = 2ir. This shows the quantization of the magnetic flux penetrating a vortex. 
The order parameter vanishes at 



•^m.n. 



i|, 



b 



for any integers m, n. The phase of (O) oc 7l(:e) rotates by 27r around each x mt „ . When 



«i 



3- 1/4 V^, 



(11.25) 



(11.26) 



the three adjoining vortices x m>n form an equilateral triangle, which is the triangular vortex lattice 
solution. Here, following [ljj, in our calculation we only focus on the case of L 2 m 2 = —1/4. Since 
different allowed values of mass only change the dimension of the condensation operator in Eq.(II.lO), 
it will not affect the Eq. (II. 17a) and also our construction of vortex solution. So we can conclude 
that for other values of m, the vortex solutions are also expected. 



III. p + ip-WAVE HOLOGRAPHIC SUPERCONDUCTOR DROPLET SOLUTION 

First we review the gravity dual theory of the non- Abclian holographic superconductor with p+i% 
wave background. The full EYM theory in 3+1 dimensional spacetime considered in Refs. (171 Hi 
has the following action 



Seym = / \f~gd A x 



M R+ i)-^ iF - Fn 



(mi) 



where #ym is the gauge coupling constant and F^ v = T a F^ v = <9 M A„ — d u A^ — i[An,A v \ is the 
field strength of the gauge field A = A^dx* = T a A°dx>*. For the SU(2) case, [T a ,T b ] = ie abc T c 
and Tr(T a T 6 ) = S ab /2, where e abc is the totally antisymmetric tensor with e 123 = 1. The Yang- 
Mills Lagrangian becomes Tr(F F P" ; ) = F^F a ^/2 with the field strength component F" = 

d^Al - d v A^ + e^A^Al. 

Working in the probe limit in which the matter fields do not backreact on the metric as in Refs. 
(171 [l8[ and taking the planar Schwarzchild-AdS ansatz, the black hole metric is the same as Eq. 
(II. 2), in which f(r) is given in Eq. (II. 3). The Hawking temperature of black hole is T = j^?-- 

Now we introduce a new coordinate z = ro/r. The metric (Eq. pi.2|l ) then becomes 

dg2 = L 2 p 2 (T) 2 + ^ + dy2) + IW 

z z z^n(z) 

where h(z) = 1 - z 3 and (3(T) = r Q /L 2 = AttT/3. 

Using the Euler-Lagrange equations, one can obtain the equations of motion for the gauge fields, 

-^dn (V^F^) + e abc A b u F c ^ = 0. (III.3) 



For the p + ip-wave backgrounds without a external magnetic field, the ansatz [17|, LL8| takes the 
following form, 

A = 4>{z, x, y)T z dt + +w(z, x, y)T x dx + w(z, x, y)T 2 dy, (III.4) 

in which we have included the spatial dependence. Here the U(l) subgroup of SU{2) generated by T 3 
is identified with the electromagnetic gauge group [13, 13 an d 4> is the electrostatic potential. Thus 
the black hole can carry charge through the condensate w, which spontaneously breaks the U(l) 
gauge symmetry below a critical temperature. This is a Higgs mechanism, but there are Goldstonc 
bosons corresponding to changing the directions of the condensate in real space or gauge space. 
They must be visible in the bulk as normal modes or (more likely) quasi-normal modes. 

In order to add a homogenous magnetic field on the boundary (where the field theory lives), we 
also need non-vanishing A 3 (z, x, y) and A^(z, x, y). Together with these non-vanishing terms above, 
the equations of motions for w, </>, A 3 A 3 are: 

2ad z (hd z w) + -(dl + d 2 Jw + ^ 2 w--w 3 ^-wd x Al + -wd y Al~-w((Al) 2 + {A 3 v ) 2 )^0 (111.5) 
a y an a a y a a ' y 

_ ad lcp - i-(9 2 + d 2 <f>) + ^w 2 = (III.6) 

an * an 



-diA% + ad z {hd z Ai) - -d x {d y A%) + -wd x w - -w 2 A y = (III.7) 

-d 2 A 3 x + ad z (hd z Al) - -d y (d x Al) - -wd x w - -w 2 A 3 x = (III.8) 

the boundary conditions are: 

(i) For w, it should be regular at the horizon. At the boundary, the asymptotic behavior of w has 
the following expression 

w=^z + -- . (111.9) 

where (O) is the condensate of the charged operator dual to the field w and is the order parameter 
for the superconductivity phase. Here we demand the constant term vanish in Eq. (III.9) since we 
require that there be no source term in the field theory action for the operator (O) [13; Oil- In fact, 
it is a requirement for the absence of such a term which in principle can be present, 
(ii) For the electromagnetic gauge fields <j>, A x and Ay, at the boundary we have 

4> = fi/P(T) -qz + --- , B(x) = d x A y - d y Al (111.10) 

in which /i is the chemical potential and q is the charge density, while -B(x) is the magnetic field. 
Obviously, at the horizon, we need = 0, and A x and A y are both regular. Now, our task is to 
solve this equation to get the information we need. 

Just as above, to exactly solve these non-linear coupled differential equations is also very difficult. 
However, as we did in the last section, we can solve the equations analytically by perturbation 
method near the upper critical magnetic field B C 2- Above B C 2 there is no condensation at any 
temperature. As in the last section, we define a deviation parameter e = (B c2 — B)/B c2 . When B is 
slightly below the upper critical magnetic field, we can expand the four fields w, (j), A x and A y as: 

w = e l/2 Wl + e 3/2 w 2 + ■ ■ ■ (III.ll) 

4> = 4> m + e<j) [l) + ■ ■ ■ (111.12) 

Al = Af) + eAf) + • • • (111.13) 

4=4(°)+e4«+--- (111.14) 

Note that all these fields are functions of x, y, z. The zeroth order solutions which correspond to the 
normal states (w = 0) with fixed chemical potential and magnetic field B c2 are 

= /i(l - z), A;f > = B c2 x, 4(°) = 0. (111.15) 

Substituting the expansion into the equation for w and making a separation of variables, we have 
w(x, y, z) = e lpv m(x 1 z;p). For a constant p, we get the equation of motion for m(x, z;p) 

2u?(l - z) 2 
[2a 2 d z (hd z ) + " \ ; }m{x,z;p) = [-d 2 x +p 2 + B 2 c2 x 2 + 3B c2 ]m(x,z;p) (111.16) 

h 

We separate the variable m as m n (x, z;p) = p n (z)j n (x;p) with a separation constant A„. The 
equations of motion for p n (z) and 7„ are: 



and 



-^(l-z) 2 +g 2 ^i) 



where X := \J2B c ix and q := /i/a are dimcnsionless. Eq. (III. 17) gives the spatial profile while 
Eq. (III. 18) gives the upper critical magnetic field. It is easy to see that the regular and bounded 
solution of Eq. (III. 17) can be expressed in terms of the Hermite function H n : 

ln (x-p) = e- x2 l A H n {X), (111.19) 

and the corresponding eigenvalue A„ is 

A„ = 2n + 4 + |-, (111.20) 

for a non-negative integer n. 

We can see that the solution of Eq. (III. 19) is independent of p, which is different from the s-wave 
model in Eq. (11.18). For the s-wave one, the spatially dependent solutions "/ n (x;p) are functions 
of p, and the vortex lattice solutions with a periodicity in x direction can be constructed by a 
superposition of different solutions for different p when n = 0. This difference leads us to conclude 
that the non-Abelian holographic superconductors cannot have vortex lattice solutions. 

These solutions are actually droplet solutions in the sense that they fall off rapidly at large |x|. 
A single droplet solution can be obtained by considering another zeroth order solution rather than 
Eq. (III. 15). We consider the following zeroth order solution: 

cj> = n(l - z), A;f > = B c2 x/2, Al^ = -B c2 y/2, (111.21) 

which satisfies the equations of motion. With this solution, after a separation of variables w(x, y, z) = 
Jn{x)"fm{y) Pm,n(z) , the solutions for the three fields are 

' ? + f)7m(*) = ¥7™(z), (IIL22) 



OX 2 8 2) " ny ' 2 

-^ + ^ + l)ln(y) = ^Jn(x), (111.23) 



_i_ {l _ z) 2 + ^M^+M] 



where X := \/2B C 2X, Y := \/2B c iy and q := fi/a are dimcnsionless. The solution for r y m {x) and 
7„(.t) are 

ln {x-p) = e- x2 '*H n {X), (111.25) 

lm (y-p) = e- Y ' 2 / & H m {X), (111.26) 
and the corresponding eigenvalue A„ is 

A m = 2m + 4, (111.27) 

A„ = 2n + 4. (111.28) 
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The order parameter for the field theory is given by the boundary value of d z w = d z p mtn j n (x)^/ m (y) 
for z = 0. p m .n is given by Eq. (III. 18), which is independent of x and y, and according to this 
equation the upper critical magnetic field has a non-zero solution satisfying the boundary conditions. 
This can be done numerically as we did in the last section, for a given q, which corresponds to a 
fixed temperature. We can also find a critical value B^/ ji 2 above which the equation has only the 
vanishing solution. The maximum upper critical magnetic field is obtained when (A„ + A m ) takes 
the minimum value (to = n = 0). The single droplet solution is also obtained when m = n = 0, 
which is Jn(x)lm(y) = e~( Y / s+x / 8 ). In Fig. 2 we give the phase diagram. 



Bc2/P 2 




VT c0 

FIG. 2. Phase diagram of the p + ip-wsve holographic superconductor in the presence of a magnetic field, 
in which we also find B C 2 oc (1 — T/T c ) at T c , just like the BCS theory. It can also be seen that the phase 
diagram is very similar to the d-wave one. 



IV. CONCLUSION AND DISCUSSION 



A d-wave and a p + ip wave holographic superconductors are studied by an analytic perturbation 
method around the upper critical magnetic field. The d-wave model has the same droplet and 
triangular vortex lattice solutions as the s-wave one, and the lattice solution is constructed by the 
superposition of droplet solutions. The p + ip-wave model has only droplet solutions, because the x 
direction property is independent of p (see Eq. (III. 19)) and so a superposition of droplet solutions 
will not give the lattice solution. According to details of our calculation, this is due to the fact 
that in the non-Abclian model the Maxwell fields appear as a U(l) subgroup of the SU(2) field and 
they do not couple with the condensed charged fields via the covariant derivative as in the s-wave 
and d-wave models. Technically, the difference between the two situations can be traced to the 
nature of the minimal coupling of the charged condensates to the gauge potential. The HTSC is 
believed to be d-wave pairing and it is indeed a type II superconductor with vortex solutions in the 
presence of a magnetic field as the holographic d-wave one. The real p + ip-wavc superconductor 
will also enter into the vortex state when a magnetic field is present, which is clearly different from 
the holographic p + ip model. When we interpret this p + zp-wave model to be a superfluidity, 
the computations described here can equally apply to this case as well with a slight modification. 
In this case, the rotation of the supcrfluid is analogous to the magnetic field. The experimental 
realization of supcrfluid droplets was given in Ref. [4J]. Although the holographic superconductor 
models have many similar behaviors as the real ones, it is still an open question to bridge the gap 
between the holographic models and the field-thcorctic method for superconductors. Finally, since 
a magnetic field can penetrate both holographic superconductors, both of them should be a type II 
superconductor. 
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